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Abstract 


Stringent limits are set on the long-lived lepton-like sector of the phenomenologi¬ 
cal minimal supersymmetric standard model (pMSSM) and the anomaly-mediated 
supersymmetry breaking (AMSB) model. The limits are derived from the results pre¬ 
sented in a recent search for long-lived charged particles in proton-proton collisions, 
based on data collected by the CMS detector at a centre-of-mass energy of 8 TeV at 
the Large Hadron Collider. In the pMSSM parameter sub-space considered, 95.9% of 
the points predicting charginos with a lifetime of at least 10 ns are excluded. These 
constraints on the pMSSM are the first obtained at the LHC. Charginos with a life¬ 
time greater than 100 ns and masses up to about 800 GeV in the AMSB model are also 
excluded. The method described can also be used to set constraints on other models. 
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1 Introduction 

We present new constraints on long-lived chargino production in the phenomenological min¬ 
imal supersymmetric standard model (pMSSM) [[Tj and on the anomaly-mediated supersym¬ 
metry breaking (AMSB) model |l2}fll. In the pMSSM, we consider the parameter sub-space for 
particle masses up to 3TeV and charginos (x^) with a mean proper decay length (ct) greater 
than 50 cm. In the AMSB model, the small mass difference between the lightest chargino and 
neutralino (^) often leads to a long chargino lifetime. 

Long-lived charged particles are predicted by various extensions of the standard model (SM) ||3- 
[Zl, such as supersymmetry (SUSY) and theories with extra dimensions dlllOl. If such par¬ 
ticles have a mass lighter than a few TeV they could be produced by the CERN LHC. The en¬ 
ergy available in the proton-proton (pp) collisions at the LHC is such that particles with mass 
>100GeV and lifetime greater than 0(1) ns could be observed with the CMS detector fTTi as 
high-momentum tracks with anomalously large rates of energy loss through ionization (d£/dx). 
These particles could also be highly penetrating such that the fraction reaching the CMS muon 
system would be sizable. The muon system could therefore be used to help in identification 
and in the measurement of the time-of-flight (TOL) of the particles. The signature described 
above is exploited in a previous CMS search [T2l, which sets the most stringent limits to date 
on a number of representative models predicting massive long-lived charged particles such as 
tau sleptons (staus), top squarks, gluinos, and leptons with an electric charge between e/3 and 
8c. 

The main thrust of this paper is to present constraints on the pMSSM and the AMSB model, 
obtained using the results from the search for heavy stable charged particles (HSCP) in proton- 
proton collisions at = 7 and 8 TeV 11121 . The method applied relies on the factorization of 
the acceptance in terms of two probabilities, calculated using the standard CMS simulation 
and reconstruction tools at 8 TeV. Moreover, the results of the acceptance calculations have 
been tabulated and made publicly available flSlI . Thus this technique, which allows the signal 
acceptance for a model with long-lived particles to be computed using the kinematic proper¬ 
ties of the particles at their production point, may in the future be used by others to evaluate 
constraints on other extensions of the SM without use of the CMS software. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter providing a field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a lead 
tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter. 
Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke of the 
magnet. The tracker measures charged particles within the pseudorapidity range 1 < 2.5 and 
provides a transverse momentum pj resolution of about 2.8% for 100 GeV particles. The ana¬ 
log readout of the tracker also enables the particle ionization energy loss to be measured with 
a resolution of about 5%. Muons are measured in the range \fi\ < 2.4, with detection planes 
based on three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. 
The muon system extends out to 11 m from the interaction point in the z direction and to 7 m 
radially. Matching tracks in the muon system to tracks measured in the silicon tracker results 
in a transverse momentum resolution between 1 and 5%, for pj values up to 1 TeV. The time 
resolution of the muon system is of the order of 1 ns. This provides a time-of-flight measure¬ 
ment that can be used to determine the inverse of the long-lived particle velocity as a fraction of 
the velocity of light (1/ j6) with a resolution of 0.065 over the full rj range I1T21 . A more detailed 
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description of the CMS detector, of the coordinate system, and of the kinematic variables can 
be found in Ref. fill . 

The first level of the CMS trigger system, composed of custom hardware processors, uses infor¬ 
mation from the calorimeters and muon detectors to select the events of interest. The high-level 
trigger processor farm further decreases the event rate from around 100 kHz to around 400 Hz 
for data storage. 

3 Estimation of signal acceptance 

Transcribing the results presented in Ref. [12l, in terms of limits on models other than those 
considered in the reference, requires the tabulation of the relevant signal acceptance. This ac¬ 
ceptance can then be used to constrain the beyond the standard model (BSM) scenarios, as 
described in the next sections. There can be significant differences in signal acceptance be¬ 
tween the models investigated in Ref. fl^ and the model to be tested. These differences arise 
because the dE/dx and TOP measurements are affected by the distribution and orientation of 
the material encountered by particles travelling within the CMS detector. The combination of 
these effects with the differences in the kinematic properties between models can result in large 
differences in signal acceptance. 

The signal acceptance can be accurately computed by fully simulating and reconstructing sig¬ 
nal events in the CMS detector and by applying the same selection criteria as adopted in 
Ref. [121 . We refer to this method of determining the signal acceptance as the "full simula¬ 
tion". This procedure requires extensive knowledge of the CMS detector and in particular use 
of CMS software that accurately models detector response to simulated signal events and em¬ 
ploys the full CMS reconstruction routines such that identical selection criteria can be used on 
simulated signal events as on data collected from collisions. 

An alternate method, which only requires information on the kinematic properties of the long- 
lived particles at their production point, is presented here. Such a method can be used if the 
efficiency for triggering and selecting events can be expressed in terms of probabilities asso¬ 
ciated with each individual long-lived particle. This is the case for models with lepton-like 
massive long-lived particles since the event selection specified in Ref. Ifl2l imposes only re¬ 
quirements on measurements performed on individual particles. The adjective "lepton-like", 
defined in Ref. 1121 , indicates particles that do not interact strongly and are therefore not subject 
to hadronization. 

The event selection requirements of Ref. Il2l are expressed in terms of measured pj, dE/dx, 
TOP, and mass values of individual particles in an event. The probability that a long-lived par¬ 
ticle in an event passes the online or offline selection requirements in Ref. I T2l can be expressed 
as a function of the true (generator-level) kinematic properties (k) of the particle: /3, rj, and pj. 
Any other equivalent set of kinematic and directional variables could also be used to express 
this dependence. The offline selection of Ref. lT2l has fixed values for pj, dE/dx, and TOP 
thresholds but uses different reconstructed mass thresholds (wJthresh) depending on the mass 
(m) of the long-lived particle in the model being tested. 

With these individual particle probabilities, the acceptance A for a model to pass the online and 
offline selections can be computed with a Monte Carlo technique by generating a large number 
of events N such that 

1 N 

^ N . X . kf). 


( 1 ) 
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where P°^{kj,k^,... ,kf^) is the probability that the event with index i containing M long- 
lived particles with true kinematic properties k],kj,... ,kf^ passes the online selection, and 
poff(j^^j^resh,kl,k?,. ■ ■ ,kf^) is the probability that the event with the same kinematic proper¬ 
ties passes the offline selection with mass threshold mthresh/ after having passed the online 
selection. Given the mass resolution of the detector, a mass threshold of wtthresh — 0.6m has 
an efficiency of about 95% for the benchmark models considered in Ref. Il2l . Consequently, 
throughout this paper, mass thresholds of 0,100, 200, and 300 GeV are used for true long-lived 
particle masses of m < 166, <330, <500, and >500 GeV, respectively. The choice of 100 GeV 
steps for the mass thresholds is made so that the information provided in Table 3 of Ref. lfT2l re¬ 
garding the background expectation and the observed count in the signal region may be used, 
thus allowing a more general application of the factorization method. In the case where only 
one long-lived particle is present in each event, the probabilities have the simplest form P°^{k) 
and thresh/^)- If each event contains two long-lived particles, the probabilities can be 

expressed using the probabilities for events with either single long-lived particle: 


P°^{k],kj) = P°^{k]) + P°^{kj) - P°^{k])P°^{kj); 

P°^^{mtbiesh,k},kj) = P°^^{mthresh,k}) + P°^^{mthresh,kj) - P°^^{mihresh,k})P°^^{mti,resh,kj)- 


( 2 ) 


The expression for events with more than two long-lived particles per event can also be ex¬ 
pressed in terms of the probabilities P°^{k) and T°^^(mthresh/^) associated with each individual 
long-lived particle. 

The probabilities P°’^{k) and T°^^(mthj.esh/^)/ derived using the full simulation as described in 
the Appendix, are evaluated using the selection requirements adopted by the "tracker+TOF" 
analysis described in Ref. fT2l . This selection, where tracks are required to be reconstructed in 
both the tracker and the muon system, has been shown to be the most sensitive to signatures 
with lepton-like long-lived particles. The probabilities thus obtained are stored in the form of 
look-up tables. 

The factorization method described above is validated by comparing the estimated signal ac¬ 
ceptance values with those obtained with full simulation for a few benchmark models predict¬ 
ing long-lived leptons. In the rest of this paper we refer to the former as "the fast technique". 
In the full simulation case, pileup due to multiple interactions per bunch crossing is also simu¬ 
lated. Agreement better than 10% between full simulation and the fast technique presented in 
this paper is generally observed for the considered values of J^thresh- 

More details on the determination of the probabilities P°^{k) and F’°^^(w^thresh/^) can be found 
in the Appendix, which also contains the details of the validation of the fast technique based 
on Eq. Q. The Appendix also explains how this technique can be used in the future to estimate 
the GMS exclusion limits for extensions of the SM other than the ones considered in this paper 
using publicly available look-up tables for P°^{k) and 

The probabilities P°^{k) and mthresh, k) are computed with simulated stable particles, but 

the method is easily extended to particles with finite lifetimes by correcting the proba¬ 

bility for the probability that the long-lived particle, with mass m, lifetime t, and momentum 
p, travels at least the distance x required to produce the minimum number of track measure¬ 
ments in the GMS muon system, as required in Ref. iT2l . The correction consists of an expo¬ 
nential factor, exp[—mx/ (rp)], to be applied to P°’^{k). This correction is not applied also to 
P°^^{mihresh,k), since it needs to be made just once for each long-lived particle. The distance v 
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only depends on the pseudorapidity of the particle: 


'9.0 m 

0.0 < |?/| < 0.8; 


10.0 m 

0.8 < \ri\ < 1.1; 

(3) 

11.0m 

\ 




These values of x ensure that the particle traverses the entire muon system before decaying. 
This choice results in a conservative estimate of the signal acceptance since it ignores the contri¬ 
bution to the acceptance from particles that decay before the end of the muon detector but still 
pass the selection. In Fig. the acceptance obtained with the fast technique is compared with 
the acceptance obtained with the full simulation of the detector, as a function of the lifetime of 
the particle. It can be seen from the ratio panels that in most cases the agreement between the 
two methods is within 10%, corresponding to the systematic uncertainty in the fast estimation 
technique used to compute the acceptance. For lifetimes less than 10 ns, the spread is somewhat 
larger, but the tendency to underestimate the acceptance is less than 15%. 




Figure 1: Signal acceptance as a function of the chargino lifetime for a benchmark model having 
a chargino of mass 100 GeV (left) and 700 GeV (right), with a mass threshold of 0 GeV and 
300 GeV, respectively. The panel below each figure shows the ratio of acceptance from the fast 
technique to the acceptance obtained from a full simulation of the detector. 

The offline event selection in Ref. IH^ includes two isolation requirements. The first is defined 
by Epj < 50 GeV, where the sum is over all tracks (except the candidate track) within a ra¬ 
dius AR = V -|- {A(p)^ = 0.3 around the candidate track. The second requirement is that 
E/p < 0.3, where E is the sum of energy deposited in the calorimeters within a radius AR = 0.3 
around the candidate track (including the candidate energy deposit) and p is the candidate 
track momentum reconstructed in the tracker. The probabilities P°^{k) and R°^^(/«thresh/^) are 
estimated with single-particle events and thus do not account for the possibility that a long- 
lived particle might fail the isolation requirements. In order to accurately model the isolation 
requirements, the following procedure, which uses generator-level information from a sim¬ 
ulation of the BSM model under test, should be used. The isolation requirements must be 
determined for each long-lived particle at the generator-level. The following conditions are 
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imposed: 

charged particles visible particles 

AR<0.3 . AR<0.3 

Pt < SOGeV and ^ E>/p < 0.3. (4) 

/ / 

The equations Q represent the generator-level equivalents of the tracker and calorimeter iso¬ 
lation requirements, respectively. The sum of the tracks pj around the long-lived candidate is 
replaced by a sum over transverse momenta of all charged particles around the direction of the 
long-lived particle. Similarly, the sum of the calorimeter energy deposits around the long-lived 
candidate is replaced by a sum over the energies of all visible particles, except the long-lived 
one, around the direction of the long-lived particle of momentum magnitude p. Muons are not 
considered in the sum over visible particles since they deposit very little energy in the calor¬ 
imeter. If a long-lived particle does not satisfy both isolation requirements expressed by Eqs. 
(j^, the corresponding factor T°^^(wJthresh/^) in Eq. (j^ is set to zero. While pileup effects are not 
taken into account in the above prescriptions, it was checked that these omissions do not have a 
significant effect on the results. The robustness of the fast estimate of the acceptance against the 
hadronic activity surrounding the long-lived particle is assessed by testing different chargino 
production mechanisms: direct pair production of charginos, pair production of gluinos that 
each decay to a heavy quark (b,t) and a chargino, pair production of gluinos that each decay to 
a light quark (u,d,s,c) and a chargino, and pair production of squarks that each decay to a quark 
and a chargino. A ~10% agreement between the fast estimate of the model acceptance and the 
full simulation predictions is observed as shown in Eig. The estimate of the acceptance for 
chargino pair production is not modified by the isolation requirement since the charginos are 
always isolated for that production mechanism. 

In the two following sections, the fast technique is used to estimate the signal acceptance of 
the CMS detector for two extensions of the SM. Given that no significant excess of events over 
the predicted backgrounds is observed in Ref. Ill2l . the signal acceptance estimated with the 
technique described above is used in calculating cross section limits of these models, at 95% 
confidence level (CE). The limits in this paper are established using the CLs approach IIT41IT51 
where p-values are computed with a hybrid Bayesian-frequentist technique IIT^ . A log-normal 
model HZl HU is adopted for the nuisance parameters. The nuisance parameters are: the ex¬ 
pected background in the signal region, the integrated luminosity, and the signal acceptance. 
The expected background and the integrated luminosity, as well as their uncertainties, are pro¬ 
vided in Ref. |12|. The uncertainty in the signal acceptance is assumed to be 25% for all signal 
models. This value results from adding 10% to the <15% uncertainty reported in Section 6 of 
Ref. mu . The additional 10% accounts for the systematic uncertainty incurred by the use of 
the fast estimation technique to compute the acceptance. Table [^summarizes the information 
needed, in addition to the signal acceptance evaluated with the fast technique, to set limits on 
a signal model predicting lepton-like charged long-lived particles. 

To demonstrate the validity of this procedure, we have used it to translate the estimate of the 
signal acceptance to the 95% GE limit on the production cross section for two of the models 
considered in Ref. [12|. The limits obtained on the pair production and inclusive production of 
staus in the context of the gauge mediated supersymmetry breaking (GMSB) model are shown 
in Table [^ for signal acceptance values estimated with full simulation and with the procedure 
described above. The limits obtained with the two techniques agree within 8%. The differences 
are due only to a small difference in the signal acceptance computed with the two different 
techniques and to the larger uncertainty assigned to the acceptance from the fast technique. 

The results in Table 1 may be compared with those shown in Table 7 and Eig. 8 of Ref. IHtI . In 
making the comparisons, it should be noted that there are important differences in the way that 



6 


3 Estimation of signai acceptance 


0 ) 

O 


80 

70 

60 

50 

40 

30 

20 

10 

0 


(8 TeV) 


= ^ ^ 1 = 

^MS Simuiation - 

— 

X* (stable, m = 200 GeV) ( 

— 

<t> 

Full simulation ^ 

L 

* 

Fast technique (w/o iso) ^ 

• 

+ 

Fast technique (with Iso) : 

L 


■ ~ 

) 


■ # 

: 


? ■ 



* 3 

- 

_^_ 

_^_1_: 


x\x. 


■ tt/bb gg^ qq 


qq 


1.2 

1 

0.8 


(8 TeV) 


80 -CMS Simulation 
70 ~- 


60 

50 

40 

30 

20 

10 

0 

1.2 

1 

0.8 


e 


X* (stable, m = 700 GeV) 

(|) Full simulation 

Fast technique (w/o Iso) 
+ Fast technique (with iso) 


x;x] 

_l_^_ L 

gg^ tt/bb gg^ qq 

qq 

- 


- 

- 


- 

- 



m 



- 


- 

- 


- 


Figure 2: Signal acceptance as a function of the chargino production mechanism for a bench¬ 
mark model having a chargino of mass 200 GeV (left) and 700 GeV (right), with a mass threshold 
of 100 GeV and 300 GeV, respectively. From left to right, the production mechanisms considered 
are: direct pair production of charginos; pair production of gluinos that each decay to a heavy 
quark (b,t) and a chargino; pair production of gluinos that each decay to a light quark (u,d,s,c) 
and a chargino; and pair production of squarks that each decay to a quark and a chargino. The 
panel below each figure shows the ratio of acceptance from the fast technique with the isolation 
requirements to the acceptance obtained from a full simulation of the detector. The estimated 
acceptance is given with and without the generator-level isolation. Pileup is present only in the 
full simulation samples. 

Table 1: Summary of the information needed to set limits on a signal model predicting lepton¬ 
like charged long-lived particles. The mass threshold, the corresponding expected background, 
and the observed numbers of events, as well as the uncertainty in the signal acceptance evalu¬ 
ated with the fast technique, are provided as a function of the long-lived particle mass. 


Mass 

(GeV) 

^thresh 

(GeV) 

Predicted 

backgrounds 

Data 

counts 

Signal Unc. 
(%) 

m < 166 

0 

44 ±9 

42 

25 

166 < m < 330 

100 

5.6 ±1.1 

7 

25 

330 < m < 500 

200 

0.56 ±0.11 

0 

25 

500 < m 

300 

0.090 ± 0.02 

0 

25 


^thresh is chosen in the two cases. In Ref. IIT^ the value was varied in steps of 10 GeV in order 
to optimise the resultant limit. In this paper a simpler approach has been followed with wJthresh 
varied in 100 GeV steps and chosen to be the largest value satisfying the condition ?«thresh < 
0.6m. The latter approach generally results in a somewhat higher estimation of the background 
and therefore a more conservative limit. In the most extreme case, the pair production of staus, 
with m = 308 GeV, the background is estimated to be 5.6 ± 1.1 events with ?«thresh = 100 GeV, 
compared with the estimate of 0.7 ± 0.1 events obtained with ?«thresh = 190 GeV in Ref. Ill2ll . 
resulting in a cross section limit that is about three times higher in the former case. Nonetheless, 
the limits agree within ~15% in almost all cases, allowing restrictive limits to be set on a general 
class of models. 
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Table 2: Signal acceptance estimated from the fast technique and with the full simulation of the 
detector, as well as the corresponding expected and observed cross section limits. Results are 
provided for both the pair production and the inclusive production of staus as predicted by the 
GMSB model. The mass threshold, the corresponding expected background and the observed 
numbers of events is also shown. 


Mass 

(GeV) 

^thresh 

(GeV) 

Predicted 

backgrounds 

Data 

counts 

Full simulation 

Acc. Exp. (fb) Ohs. (fb) 

Fast technique 

Acc. Exp. (fb) Ohs. (fb) 

Pair production of staus 

126 

0 

44 ±9 

42 

0.24 

4.38 

4.11 

0.24 

4.53 

4.24 

156 

0 

44 ±9 

42 

0.28 

3.66 

3.43 

0.29 

3.81 

3.57 

200 

100 

5.6 ±1.1 

7 

0.34 

1.06 

1.28 

0.35 

1.08 

1.30 

247 

100 

5.6 ±1.1 

7 

0.40 

0.90 

1.09 

0.40 

0.93 

1.13 

308 

100 

5.6 ±1.1 

7 

0.46 

0.77 

0.93 

0.47 

0.79 

0.96 

370 

200 

0.56 ±0.11 

0 

0.53 

0.41 

0.31 

0.53 

0.42 

0.32 

494 

200 

0.56 ±0.11 

0 

0.61 

0.36 

0.27 

0.62 

0.37 

0.28 

745 

300 

0.09 ± 0.02 

0 

0.66 

0.24 

0.24 

0.67 

0.25 

0.24 

1029 

300 

0.09 ± 0.02 

0 

0.58 

0.28 

0.27 

0.59 

0.28 

0.27 

Inclusive production of staus 

126 

0 

44 ±9 

42 

0.25 

4.22 

3.95 

0.25 

4.43 

4.15 

156 

0 

44 ±9 

42 

0.32 

3.21 

3.01 

0.32 

3.38 

3.16 

200 

100 

5.6 ±1.1 

7 

0.41 

0.87 

1.05 

0.42 

0.90 

1.09 

247 

100 

5.6 ±1.1 

7 

0.50 

0.72 

0.87 

0.50 

0.76 

0.91 

308 

100 

5.6 ±1.1 

7 

0.56 

0.64 

0.77 

0.56 

0.67 

0.81 

370 

200 

0.56 ±0.11 

0 

0.60 

0.36 

0.27 

0.60 

0.37 

0.28 

494 

200 

0.56 ±0.11 

0 

0.66 

0.33 

0.25 

0.65 

0.35 

0.26 

745 

300 

0.09 ± 0.02 

0 

0.67 

0.24 

0.23 

0.67 

0.25 

0.24 

1029 

300 

0.09 ± 0.02 

0 

0.58 

0.28 

0.27 

0.58 

0.29 

0.28 


Because the online selection in Ref. IH^ uses a missing transverse energy trigger in com¬ 

bination with a single-muon trigger, there is one caveat to the proposed factorization method: 
the efficiency of the trigger cannot be modeled accurately in terms of single long-lived 
particle kinematic properties. Accounting for the presence of other undetectable particles us¬ 
ing a Monte Carlo method would not help because £™^®® often depends significantly on detector 
effects due to the other particles. The assumption that the £Jp‘®® trigger adds negligibly to the 
event selection performed by the muon trigger must therefore be satisfied in order to apply the 
method to a given signal. Deviations from this assumption would result in an underestimation 
of the signal acceptance. The assumption is satisfied by models with lepton-like long-lived par¬ 
ticles. Models with long-lived colored particles, such as top squarks or gluinos, do not satisfy 
this condition and thus cannot currently be tested with the technique presented in this paper. 
Long-lived colored particles hadronize in color singlet bound states fl^ that could interact 
with the detector material leading to complex situations, detailed in Refs, [ill |20l, and induc¬ 
ing significant instrumental £™‘®®. For instance, pair-produced colored long-lived particles may 
hadronize to a charged and a neutral hadron. In this case the £™'®® is strongly modified by the 
presence of the neutral hadron since it is not visible in the tracker and deposits only 0(1) GeV 
in the calorimeter. Moreover, the interactions of the color singlet bound states with matter may 
lead to a modification of the electric charge of the bound states 1191 . In this case, the hadron 
containing a long-lived colored particle can be electrically charged at production, but neutral 
in the muon system and therefore fail the muon reconstruction. The probabilities P°^{k) and 
thresh/^) do not account for the possible modification of the electric charge experienced 
by the hadron-like particles. 
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4 Constraints on the pMSSM 

The pMSSM is a 19-parameter realization of the minimal supersymmetric model (MSSM) IT] 
that captures most of the phenomenological features of the R-parity conserving MSSM. The 
free parameters of the pMSSM, in addition to the SM parameters, are: (1) the gaugino mass 
parameters Mi, M 2 , and M 3 ; (2) the ratio of the vacuum expectation values of the two Higgs 
doublets tan ji = V 2 /v\) (3) the higgsino mass parameter ^ and the pseudoscalar Higgs mass 
mfi,) (4) the 10 sfermion mass parameters m^, where F = Qi,lJi,Di,Li,£i,Q 3 , 113 , 03 , 1 . 3 , £3 
(imposing degeneracy of the first two generations niQ^ = ...); and 5) the 

trilinear couplings At, A^, and Aj. To minimize the theoretical uncertainties in the Higgs sector, 
these parameters are defined at the electroweak symmetry breaking scale. 

In the pMSSM, all MSSM parameters are allowed to vary freely, subject to the requirement 
that the model is consistent with some basic constraints: first, the sparticle spectrum must be 
free of tachyons and cannot lead to color or charge breaking minima in the scalar potential. 
We also require that the model is consistent with electroweak symmetry breaking and that 
the Higgs potential is bounded from below. Finally, in this study, we also require the lightest 
SUSY particle to be the lightest neutralino (^). Furthermore, for practical reasons, we limit 
our study to the pMSSM sub-space chosen to cover sparticle masses up to about 3 TeV. Table 
summarizes the boundaries of the considered sub-space. 

Table 3: The pMSSM parameter space used in the scan. 


-3 < 

Ml, M 2 

< 3 TeV 

0 < 

M 3 

< 3 TeV 

-3 < 

F 

< 3 TeV 

0 < 

niA 

< 3 TeV 

2 < 

tan/3 

<60 

0 < 

Qi, 2 / Qa 

< 3 TeV 

0 < 

Hi,2,U3 

< 3 TeV 

0 < 

F> 1 , 2 / C >3 

< 3 TeV 

0 < 

kl, 2 / ^3 

< 3 TeV 

0 < 

^1,2/ £3 

< 3 TeV 

-7 < 

At, A\,, Ax 

< 7TeV 


In Ref. I2T]| , some previous CMS published results were reinterpreted in the context of this 
pMSSM sub-space. This analysis, however, did not consider the region of parameter space in 
which long-lived charginos are predicted. Using the technique described in Section we can 
extend the results of Ref. Kill to regions of the parameter space leading to long-lived particles. 

We sample 20 million points in a pMSSM parameter space from a prior probability density 
function that encodes results from indirect SUSY searches and pre-LHC searches as done in 
Ref. |21]. From this set we select 7205 points with a Higgs boson mass in the range 120 < 
m-^ < 130 GeV and predicting long-lived (ct > 50 cm) charginos. Tightening further the mass 
window to 123 < m-^ <128 GeV to reflect more recent constraints from the Higgs boson mass 
measurements 122112411 reduces the size of the subset by 45% but does not further constrain the 
chargino mass or lifetime. We therefore use the 120 < < 130 GeV window in order to 

minimize the statistical uncertainty. 

For each of the 7205 points in this subset, we have generated 10000 events using PYTHIA 
v6.426 l25l . Both direct pair production of charginos and indirect chargino production through 
the decay of heavier SUSY particles were considered for this study. The generated events have 
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been used to evaluate the signal acceptance of the HSCP search, given the chargino kinematic 
properties predicted by PYTHIA for the considered pMSSM sub-space. 

The fast technique is used to obtain acceptance values expected to be in good agreement with 
the full simulation prediction. The predicted signal acceptance is then used to compute 95% CL 
limits on the 7205 analysed pMSSM parameter points. A parameter point is excluded if the ob¬ 
served limit obtained on the cross section is less than the theoretical prediction at leading-order 
as calculated by PYTHIA. The use of leading-order instead of next-to-leading-order theoretical 
cross section is driven by practical considerations given the large number of parameter points 
considered. 


Figure shows the fraction of parameter points excluded as a function of the chargino life¬ 
time. The fraction of excluded model points with a chargino lifetime longer than 1000 ns 
(10 ns) is 100.0% (95.9%). Although these values depend on the random point sampling in 
the 19-dimensional pMSSM parameter space, it is remarkable that a high fraction of the points 
predicting long-lived charginos are excluded. 

Figure]^ shows the number of parameter points predicted and excluded by the analysis of the 
results obtained in Ref. |[12H as a function of the chargino mass, chargino lifetime, and the mass 
difference between the chargino and the neutralino. 




Figure 3: (left) Number of pMSSM points, in the sub-space covering sparticle masses up to 
about 3 TeV, that are excluded at a 95% CL (hatched red) or allowed (white) as a function of 
the chargino lifetime, (right) Enlargement of the long-lived region. The bottom panel shows 
the fraction of pMSSM points excluded by the analysis based on the results from the HSCF 
search Ill2l . 


5 Constraints on the AMSB model 

In the AMSB model 121-01 the lightest chargino and neutralino are almost degenerate (Piy± — 
< 1 GeV), where the neutralino is the lightest SUSY particle. In this model, the chargino 
lifetime, expected to be of the order of a nanosecond or larger, is determined by the mass 
splitting with the neutralino. 
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Figure 4: Number of pMSSM parameter points in the sub-space covering sparticle masses up to 
about 3 TeV shown as a function of the chargino mass and (upper row) of the mass difference 
between the chargino and the neutralino, and (lower row) chargino lifetime. The left panels 
show the entire set of points considered while the right panels show the set of points excluded 
by the analysis based on the results from the HSCP search Ill2l . 
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Previous searches for AMSB charginos Il2^l27ll looked for a chargino decaying within the track¬ 
ing volume into a neutralino and a soft charged pion. The pion has a momentum of ~ 100 
MeV and is generally not reconstructed. The experimental chargino signature therefore takes 
the form of a disappearing track inside the tracking system. The main limitation of that search 
technique is that the sensitivity drops quickly as the lifetime increases because the chargino is 
required to decay within the tracking region. 

In contrast, the sensitivity of the search for HSCP lfT2ll is maximal when the chargino decay 
occurs outside the detector. These two searches are therefore complementary. In this context, 
the fast technique discussed in the previous sections can be used to assess the limits set by 
Ref. 1121 on long-lived charginos in the AMSB model. 

The minimal version of the AMSB model is fully characterized by four parameters: the ratio 
of Higgs doublet vacuum expectation values at the electroweak scale, the sign of the higgsino 
mass term, the universal scalar mass (mo), and the gravitino mass {m^/ 2 ), which dictates the 
value of the chargino mass. The values of the first two parameters are set to tan /3 = 5 and 
}i > 0. The scalar mass is set to a large value (1 TeV) in order to prevent the appearance of 
tachyonic sleptons. Gravitino masses ranging from 3.5 to 32 TeV are used in order to scan 
chargino masses from 100 to 900 GeV. 

Samples of simulated charginos with masses from 100 to 900 GeV and lifetimes from 1 ns to 
10 jis are produced with pythia v6.426. The SUSY mass spectrum and the decay tables are 
calculated using ISASUSY with ISAJET v7.80 |28l. For each sample, 10000 events are generated 
using an inclusive SUSY production and the acceptance of the search for long-lived particles 
is estimated using the technique described in Section The estimated signal acceptance is 
shown in Fig. (left). The acceptance is reduced for short lifetimes because the probability that 
a particle reaches the muon system before decaying, is exponentially smaller. As explained in 
Section]^ a systematic uncertainty of 25% is assigned to the signal acceptance. A point in the 
mass-lifetime parameter space of the AMSB model is considered to be excluded when the 95% 
GL observed limit on the cross section is lower than the leading-order theoretical cross section. 
The excluded region in this plane is shown in Fig. (right). These results extend those from 
previous searches at LHG experiments ||2^|27| by excluding charginos with lifetimes >100 ns 
up to masses of about 800 GeV. While the signal acceptance remains nearly constant over a 
wide mass range, heavier charginos cannot be excluded because of their smaller production 
cross section. The sensitivity of the search for HSGF lll^ is limited to charginos with an average 
lifetime > 3 ns while previous searches based on short track signatures are sensitive to lifetimes 
of ~ 0.1ns liilzzi. 

6 Summary 

The results of the GMS search for long-lived charged particles have been analysed to set con¬ 
straints on the phenomenological minimal supersymmetric standard model (pMSSM) and on 
the anomaly-mediated SUSY breaking model (AMSB), both of which predict the existence of 
long-lived massive particles in certain regions of their parameter space. A novel technique 
for estimating the signal acceptance with an accuracy of 10% is presented. This technique 
only uses generator-level information, while the integrated luminosity, the expected standard 
model background yields, and the corresponding uncertainties are taken from a previous GMS 
search Iil2l . The technique and the tabulated probabilities, available as supplementary material 
to this paper (131, can be used by others to estimate the GMS exclusion limits for different mod¬ 
els predicting long-lived lepton-like particles. In the context of the AMSB model, charginos 
with lifetimes > 100 ns (3 ns) and masses up to about 800 GeV (100 GeV) are excluded at 95% 
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Figure 5: (left) Signal acceptance as a function of chargino mass for the AMSB model as pre¬ 
dicted by the fast technique, (right) Observed and expected excluded region on the chargino 
mass and lifetime parameter space in the context of the AMSB model with tan j6 = 5 and > 0. 
The excluded region is indicated by the hatched area. 

confidence level. The most stringent limits to date are set on the long-lived sector of the pMSSM 
sub-space that covers SUSY particle masses up to about 3 TeV. In this sub-space, 95.9% (100%) 
of the points with a chargino lifetime r > 10 ns (1000 ns) are excluded by the present analysis 
of the results from the CMS search in Ref. [13. These are the first constraints on the pMSSM 
obtained at the LHC. 
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A Details of the fast technique 

The probabilities P°^{k) and P°^^(mthresh/^)/ introduced in Section]^ are computed using sam¬ 
ples of single long-lived particles, uniformly distributed in y and j6, and produced in PYTHIA 
v6.426 1251 . Twenty samples, each containing one million stau particles, were produced for the 
following long-lived particle masses: 100, 126, 156, 200, 247, 308, 370, 432, 494, 557, 651, 745, 
871, 1029, 1200, 1400, 1600, 1800, 2000, and 2500 GeV. These generated events were then pro¬ 
cessed by the full GMS simulation and reconstruction software. The full simulation includes 
pileup effects. The selection requirements adopted by the "tracker+TOF" analysis described in 
Ref. were then applied to the reconstructed events in order to evaluate the probability that 
a particle with kinematics k can pass the selection criteria. 


The probabilities P°^{k) and thresh/^) were evaluated in 3D bins of \fi\, /3, and pj. For 

\i]\ we consider the following bin boundaries: 0.0, 0.25, 0.50, 0.75, 1.00, 1.10, 1.125, 1.50, 1.75, 
2.00, 2.10. The granularity is finer around the transition between the tracker barrel and the 
tracker endcap to better model the efficiency in this region. Since the detector is symmetric in 
rj, the probabilities for the positive and negative y regions were averaged in order to reduce 
the statistical uncertainties. A constant bin width of 0.05 is considered between 0.0 and 1.0 for 
the binning in j6. For the Firming along pj the following bin boundaries are considered 5, 50, 


60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 220, 240, 260, 280, 300, 325, 
350,375,400,425,450,475,500,550,600,650, 700, 750,800,850,900,950,1000,1100,1200,1300, 
1400, 1500, 1600, 1700, 1800, 1900, 2000, 2250, 2500, ooGeV. Figures \kj\ and \A^ show 
the estimated values of P°’^{k), P°^^(0GeV,k), and P°^^(300GeV,k), respectively. The complete 
evaluation of P°^{k) and 4*°^^(mthresh/^) for the bin ranges previously described is available as 
supplemental material to this paper 


The fast technique is validated by comparing the estimated signal acceptance values with those 
obtained from the full simulation for three benchmark models predicting long-lived leptons: 
pair production of staus, inclusive production of staus in the context of a gauge mediated sym¬ 
metry breaking (GMSB) model, and pair production of long-lived leptons with a null weak 
isospin |[29| . These three models, having significantly different kinematic properties, are stud¬ 
ied in Ref. [121. 


The signal acceptance obtained with the two methods is shown in Figs. A.4 and A.5 for the 
benchmark models considered for two thresholds on the reconstructed mass, wJthresh- The ac¬ 
ceptance computed with the full simulation and with the technique presented in this paper 
agrees within 10%. The agreement between the two techniques is worse when the mass thresh¬ 
old is close to the true mass of the long-lived particle. The larger disagreement close to the 
mass threshold is a consequence of the coarse binning of probability F’°^^(?«thresh/ k), which only 
partially reflects the sharp variation in probability with the reconstructed values of pj, y, and 
mass. The reconstructed particle mass is generally lower than the actual particle mass because 
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A Details of the fast technique 


of a cutoff applied in the data acquisition electronics of the tracker detector. In Ref. il2l , the 
threshold on the reconstructed mass is optimized for each signal model considered. In this 
paper, for the sake of simplicity, the threshold on the reconstructed mass used in the fast tech¬ 


nique is set at 60% of the true particle mass. The hatched area in Figs. A.4 and A.5 indicates 
the range not satisfying the requirement on the reconstructed mass threshold. Outside of this 
region the estimated signal acceptance is always compatible within 10% with the results from 
full simulation. 


Several points of the pMSSM subset have also been fully simulated, as described in Ref. 
in order to further validate the fast estimation of the acceptance in the context of pMSSM. 
The acceptances estimated using the full simulation prediction and the proposed technique are 


compared in Fig. A.6 for points leading to quasi-stable (ct >100 m) charginos. An agreement 
at the level of 10% is observed. 
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Figure A.l: Values taken by the probability as a function of the true particle-variables 

Pj, /3, and \r]\. The binning in pj in these figures is coarser than the one used in the study 
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Figure A.2: Values taken by the probability P°^^{0,k) as a function of the true particle-variables 
Pj, j6, and \rj\. The binning in pj in these figures is coarser than the one used in the study 


(r^hresh = 0 GeV, k) = 0 GeV, k) = 0 GeV, k) 
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Figure A.3: Values taken by the probability P°^^(300GeV,A:) as a function of the true particle- 
variables Pj, j6 , and \rj\. The binning in px in these figures is coarser than the one used in the 
study 


Mhresh = 300 GeV, k) = 300 GeV, k) = 300 GeV, k) 
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Figure A.4: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right). The upper 
and lower sets of distributions show the pair production and inclusive production of staus as 
predicted by the GMSB model, respectively. The panel below each figure shows the ratio of 
acceptance from the fast technique to the acceptance obtained from a full simulation of the de¬ 
tector. The hatched area indicates the range not satisfying the requirement on the reconstructed 
mass threshold. 
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Figure A.5: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right). The tested 
model is the pair production of leptons with no weak isospin. The panel below each figure 
shows the ratio of acceptance from the fast technique to the acceptance obtained from a full 
simulation of the detector. The hatched area indicates the range not satisfying the requirement 
on the reconstructed mass threshold. 
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Figure A.6: Signal acceptance for a mass threshold of 0 (left) and 300 GeV (right) on a few 
representative pMSSM points predicting quasi-stable (ct > 100 m) charginos. The panel below 
each figure shows the ratio of acceptance from the fast technique to the acceptance obtained 
from a full simulation of the detector. The hatched area indicates the range not satisfying the 
requirement on the reconstructed mass threshold. 
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